WORLD INTELLECTUAL PROPERTY ORGANIZATION 

International Bureau 




PCT 

INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 



(51) International Patent 

G01V 3/20, 3/28, 3/38 



Al 



(11) International Publication Number: 



WO 95/03557 



(43) International Publication Date: 2 February 1995 (02.02.95) 



(21) International Application Number: 

(22) International Filing Date: 



PCT/US94/07991 
19 My 1994 (19.07.94) 



(30) Priority Data: 
08/095394 



21 July 1993 (21.07.93) 



US 



(71) Applicant: WESTERN ATLAS INTERNATIONAL, INC. 

[US/US]; 10205 Westheimcr Road, Houston, TX 77042- 
4299 (US). 

(72) Inventors: STRACK, Kurt, M; 10811 Holly Springs, Houston, 

TX 77042 (US). TABAROVSKY, Leonty, A; 13739 Piping 
Rock, Houston, TX 77077 (US). PAYTON, Christopher, C; 
2424 S. Voss #311, Houston, TX 77057 (US). 

(74) Agent: SPRINGS, Danyl, M; Western Adas International, 
Inc., 10205 Westheimer Road, Houston, TX 77042-4299 
(US). 



(81) Designated States: CA, NO, European patent (AT, BE, CH, 
DE, DK, ES, FR, GB, GR, IE, IT, LU, MC, NL, PT, SE). 



Published 

With international search report 

Before the expiration of the time limit for amending the 
claims and to be republished in the event of the receipt of 
amendments. 



(54) Title: METHOD OF DETERMINING FORMATION RESISTIVITY UTILIZING COMBINED MEASUREMENTS OF INDUCTIVE 
AND GALVANIC LOGGING INSTRUMENTS 



138 




138 



140 



(57) Abstract 

A method of determining selected parameters of an earth formation (132) surrounding a borehole (135) by first obtaining at least one 
induction logging measurements of the selected parameters in a first predetermined volume (138-138*) of the formation (132) surrounding 
the borehole (135) having known first radial and vertical dimensions, then obtaining at least one galvanic logging measurement of the 
identical selected parameters in a second predetermined volume (136) of the formation (132) surrounding the borehole (135) having known 
second radial and vertical dimensions that overlap the first radial and vertical dimensions of the first predetermined volume, whereby the 
overlapping volumes form a representative common volume (140) of the formation (132), and then combining the induction and galvanic 
logging measurements using an inversion technique to obtain a measurement of the selected parameters of the earth formation (132) 
surrounding the borehole (135) in the representative common volume (140) of the formation (132). 
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METHOD OF DETERMINING FORMATION RESISTIVITY UTILIZING 
COMBINED MEASUREMENTS OF INDUCTIVE AND GALVANIC 

LOGGING INSTRUMENTS 

This invention relates to methods of determining earth formation parameters and 
more particularly relates to a method of measuring formation resistivity surrounding 
a borehole utilizing combined measurements of inductive and galvanic logging 
instruments. 

Modern conventional electrical resistivity measurements are grouped into two 
classes, those injecting electrical currents into the formation by means of electrodes 
(galvanic logging devices, including lateral, spherically focused, and normal devices) 
and those using coils (induction logging devices) for creating eddy currents in the 
formation. The galvanic logging methods are really just developments of the original 
electrode instrument methods invented by the Schlumberger brothers in the 1920' s. (L. 
Allaud and M. Martin, "Schlumberger, the History of a Technique," Wiley, New 
York, 1977). The induction logging methods and devices were created by Henri- 
Georges Doll in the 1940's. (H.G. Doll, Pet. Trans. AIME. 186, 148, 1947). 

For the induction logging device the signal measured from a particular 
formation zone is inversely related to the resistivities in the formation around the 
borehole and to the resistivities within the borehole. For the galvanic logging devices 
the signal measured is non-linearly related to the resistivities and the resistivity 
contrasts of the borehole and the formations surrounding the borehole. At the 
borehole/formation/invaded zone boundaries, electrical charges arise for galvanic 
logging devices, while for induction logging devices, the effect is due to the induced 
current in the borehole/formation/invaded zone media. In the field of electromagnetic 
surveying, the different situations are usually described as different modes, the 
transverse electric (TE) mode in the case of the induction devices, and the galvanic or 
transverse magnetic (TM) mode for the galvanic (electrode) devices. 
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In the field of surface electromagnetic surveying of the earth, the joint 
measurement and interpretation of inductive and galvanic (TE and TM) techniques has 
been investigated in recent years. (T. Eadie, "Detection of Hydrocarbon 
Accumulations by Surface Electrical Methods - A Feasibility Study," Master's Thesis, 
University of Toronto, 1980). Based on Eadie's work and the initial application of a 
numerical interpretation technique referred to as "joint inversion" (K. Vozoff and 
D.L.B. Jupp, "Joint Inversion of Geophysical Data," Geophvs. J.R. astr. Soc . v. 42, 
pp. 977-991, 1975; D.L.B. Jupp and K. Vozoff, "Resolving Anisotrophy in Layered 
media by Joint Inversion," Geophvs. Pros. v. 25, pp. 460-470, 1977; and A.P. Raiche, 
D.L.B. Jupp, H. Rutters and K. Vozoff, "The Joint Use of Coincident Loop Transient 
Electromagnetic and Schlumberger Sounding to Resolve Layered Structures," 
Geophysics , v. 50, n. 10, pp. 1618-1627, 1985) led to the application of the combined 
modes (inductive and galvanic) for surface geophysical prospecting (K.-M. Strack, T. 
Hansteih, K. LeBrocq, D.C. Moss, K. Vozoff and P.A. Wolfgram, "Case Histories 
of LOTEM Surveys in Hydrocarbon Prospective Areas," First Break, v. 7, n. 12, pp. 
467-476, 1989; K.-M. Strack, "Exploration With Deep Transient Electromagnetics," 
Elsevier, 1992, 373 pp.). However, for borehole applications the combination of the 
two modes (inductive and galvanic) to characterize the same volume of formation has 
not been developed. Recently, a paper has indicated that the inversion of induction 
tool responses alone has become available as a wellsite product as an improved 
interpretation method. (A.Q. Howard, "A New Invasion Model for Resistivity Log 
Interpretation," The Log Analyst , pp. 97-110, March-April 1992). Historically, 
induction and galvanic logging measurements investigating the same part of the 
formation are usually used in different borehole environments (mud conductivity) 
(Ellis, "Well Logging for Earth Scientists," Elsevier, (1987)) and for different relative 
formation resistivities (D.E. Johnson and K.E. Pile, "Well Logging For the 
Nontechnical Person," pp. 68-70, PennWell, 1988). Based on these considerations, 
either the induction logging or galvanic logging device is chosen for the 
borehole/formation resistivity conditions. Rarely are both logging devices (induction 
and galvanic) which investigate the same or similar volumes of formation run together 
in the borehole because of the excessive tool string length required and the increased 
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cost. Although sometimes induction and galvanic devices are run in the same borehole, 
usually in separate runs, their measurements are not generally used to produce a 
combined interpreted response because either one of the logs produced by the two 
devices is deemed to exhibit the most applicable response characteristics to each 
particular section of the borehole and the formation. 

The method utilizes the combined measurements from both induction and 
galvanic logging instruments to obtain an improved definition of the electrical 
parameters of the formation surrounding a borehole. 

The measurement data sets can be derived by using borehole logging devices 
giving induction and galvanic logging type responses which could have the following 
nature: 

- two or more separate tools used in either separate logging runs or as a 
combined tool string. 

- one combined tool measuring the inductive and galvanic response of the 
medium in parallel (simultaneously or separated by a short interval). 

- The combination can be performed by (i) one mandrel, (ii) two or more 
mandrels on one tool string, or (iii) separate tools run separately in the borehole 
even from different tool families as long as they characterize the induction and 
galvanic response of the medium. 

The same mathematical model for both types of data sets is used in iterative 
optimization techniques which yield one optimized model consistent with both data sets 
for a given situation. The uncertainty in the definition of the model output parameters 
is improved through the combined usage of the two (or more) different data sets, and 
thus the precision and accuracy of the resulting formation resistivity value is improved 
(which is the key indicator of the interpretation quality.) 

The optimization can be done in several ways including using iterative 
numerical techniques, or graphical techniques including correction charts, or even 
direct imaging techniques. The convergence to the final model is data and model 
parameters driven. The reliability of the improvement can be evaluated by using error 
bounds or confidence levels or sensitivities to the model parameters. 
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In order that the manner in which the above-recited principles and features of 
the invention are attained can be understood in detail, a more particular description of 
the invention may be had by reference to specific embodiments thereof which are 
illustrated in the accompanying drawings, which drawings form a part of this 
specification. 

In the drawings: 

Fig. 1 is a simplified schematic drawing depicting a two-coil induction logging 
instrument including a diagrammatic presentation of the measurement technique of the 
induced electromagnetic field introduced into a formation. 

Fig. 2 is a simplified schematic drawing of a typical three-electrode lateral 
logging instrument including a diagrammatic presentation of the measurement technique 
of the electrical current introduced into a formation. 

Fig. 3 is a simplified schematic drawing of a typical seven-electrode dual lateral 
logging instrument including a diagrammatic presentation of the measurement technique 
of the electrical current introduced into a formation. 

Fig. 4 is a diagram illustrating the typical environmental and formation 
conditions for currently selecting the use of a lateral or induction instrument when 
logging. 

Fig. 5 is a stylized schematic diagram of an assumed earth model formation 
penetrated by a borehole and showing various formation medium related parameters. 

Fig. 6 is an example of a log made using a combination of a selected dual phase 
induction logging (DPIL) instrument and a selected shallow focused laterolog (SFL) 
instrument showing the DPIL deep and medium responses and the SFL shallow 
response. 

Fig. 7 is a chart for providing a method of determining Rt from the DPIL log 

readings where ^ < R^. 

Fig. 8 is a chart for providing a method of determining Rt from the DPIL log 

readings where R t >R x0 . 

Fig. 9 is an example of a log made using a selected dual laterolog (DLL) 

instrument showing the DLL deep and shallow responses. 

Fig. 10 is a chart for providing a method of determining Rt from DLL log 
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readings as shown in Fig. 9. 

Fig. 11 is an example of a log made using a combination of a selected dual 
phase induction logging (DPIL) instrument and a selected micro-laterolog (MLL) 
instrument showing the DLL deep and SLL shallow responses and the MLL response. 

Fig. 12 is a schematic drawing of a selected combination logging tool including 
both a selected high vertical resolution lateral logging instrument for measuring thin- 
bed resistivity (TBRT) and a micro-laterolog (MLL) instrument. 

Fig. 13 is an example of a log made using a selected high vertical resolution 
lateral logging instrument for measuring thin-bed resistivity (TBRT). 

Fig. 14 is a diagram showing two sets of a plurality of response curves for each 
of the selected TBRT and MLL instruments shown in Fig. 12 for different values of 

the ratio of R t to R^. 

Fig. 15 is a stylized schematic diagram of a borehole penetrating on earth 
formation showing a common volume of formation surrounding the borehole that is 
investigated by selected induction and galvanic logging instruments. 

Fig. 16 is a stylized schematic diagram of a borehole penetrating on earth 
formation showing common multiple volumes of formation surrounding the borehole 
at radially spaced distances that are investigated by selected induction and galvanic 
logging instruments. 

Fig. 17 is a stylized schematic diagram of a borehole penetrating on earth 
formation showing common multiple volumes of formation surrounding the borehole 
at vertically spaced distances that are investigated by selected induction and galvanic 

logging instruments. 

Fig. 18 is a diagram showing a three-dimensional representation of an 
illustrative conceptual example of ellipsoids of confidence for percentage deviations of 
the parameters Rt, Rxo, and Di for the earth model shown in Fig. 5 projected to the 
axes of the earth model parameters for the induction and lateral logging instruments 

referred to in Figs. 6-14. 

Fig. 19 is a diagram showing the ellipsoids of confidence for percentage 
deviations of the earth model parameters R^ R xo , and D { for a first selected earth model 
using preselected combinations of the induction and lateral logging instruments referred 
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Fig. 20 is a diagram showing a plot of the formation resistivity vs. the radial 
distance or depth into the formation for reflecting the confidence level envelopes of the 
ellipsoids of confidence shown in Fig. 19 for the first selected earth model using the 
preselected combinations of the induction and lateral logging instruments. 

Fig. 21 is a diagram showing the ellipsoids of confidence for percentage 
deviations of the earth model parameters Rt, R x0 , and D t for a second selected earth 
model using preselected combinations of the induction and lateral logging instruments. 

Fig. 22 is a diagram showing a plot of the formation resistivity vs. the radial 
distance or depth into the formation for reflecting the confidence level envelopes of the 
ellipsoids of confidence shown in Fig. 21 for the second selected earth model using the 
preselected combinations of the induction and lateral logging instruments. 

Fig. 23 is a diagram showing the ellipsoids of confidence for percentage 
deviations of the earth model parameters Rt» R xo , and Dj for a third selected earth 
model using preselected combinations of the induction and lateral logging instruments. 

Fig. 24 is a diagram showing a plot of the formation resistivity vs. the radial 
distance or depth into the formation for reflecting the confidence level envelopes of the 
ellipsoids of confidence shown in Fig. 23 for the third selected earth model using the 
preselected combinations of the induction and lateral logging instruments. 

Fig. 25 is a diagram showing the ellipsoids of confidence for percentage 
deviations of the earth model parameters R,, R x0 , and D, for a fourth selected earth 
model using preselected combinations of the induction and lateral logging instruments. 

Fig. 26 is a diagram showing a plot of the formation resistivity vs. the radial 
distance or depth into the formation for reflecting the confidence level envelopes of the 
ellipsoids of confidence shown in Fig. 25 for the fourth selected earth model using the 
preselected combinations of the induction and lateral logging instruments. 

Fig. 27 is a functional block diagram schematic of the steps of an example 
inverse modelling technique used in combining field induction and galvanic logging 
data to determine a desired formation resistivity parameter in accordance with this 
invention. 
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Referring now to Figs. 1-3, a brief description of the construction and operation 
of typical simplified induction and galvanic electrical resistivity logging instruments 
will be given. In Fig. 1, a borehole 30 is shown penetrating an earth formation 32 and 
in which is suspended an elongated induction logging instrument 34 by means of a 
wireline cable 36 communicating with the earth's surface (not shown). The logging 
tool 34 is raised and lowered by means of the cable 36 using conventional winch-type 
hoist trucks/skids (not shown), and command and data signals are transmitted to and 
from the tool to the surface through the cable. Spaced longitudinally in the tool 
housing 34 is a transmitter coil 40 interconnected to a transmitter oscillator 42 for 
generating a transmitter signal of a desired frequency and a receiver coil 44 
interconnected to a receiver circuit 46. The transmitter coil 40 is energized by the 
transmitter oscillator 42 which sets up a magnetic field around the tool and extending 
into the formation. The constantly changing magnetic field causes currents to flow 
through the formation concentric to the axis of the tool 34. The currents, called 
ground loops, are shown diagrammaticaily as a ring-shaped section 48 of the formation 
32 and are proportional to the conductivity of the formation. The ground loops 
currents at 48 alternate at the same frequency as the magnetic field and the transmitter 
current flowing through the transmitter coil 40. The ground loop currents 48 also set 
up secondary magnetic fields of their own, which interact with the receiver coil 44 to 
cause an electrical current to flow in the receiver coil 44. The amount of the current 
flowing in the receiver coil 44 is related to the ground loop 48 currents and therefore 
to the conductivity of the formation. The currents in the receiver coil 44 are detected 
by the receiver circuit 46 and transmitted to the surface for processing and recording 
(not shown) as either a conductivity measurement C or a resistivity measurement R, 
where R = 1/C. 

Fig. 2 shows a simple three-electrode "focused" electrical galvanic resistivity 
instrument 56 of the type known as a lateral logging device disposed in a borehole 30 
penetrating an earth formation 32 and suspended therein by a wireline cable 58. The 
logging tool 56 is raised and lowered by means of the cable 58 by conventional means 
as previously described (not shown) , and command and data signals are transmitted to 
and from the tool to the surface through the cable 58. The tool 56 includes an upper 
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"guard" or "focusing" electrode 60 (Al), a lower "guard" or "focusing" electrode 62 
(A2) and a "center" or "measuring" electrode 64 (AO). A constant current is emitted 
from the center or measuring electrode 64 (AO). The electrical potential of the pair of 
guard or focusing electrodes 60 (Al) and 62 (A2) are maintained equal in order that 
the current emitted from electrode 64 (AO) is forced out "laterally" into the formation 
and is "focused" in the section shown by the shaded area 66. 

Another more complex "focused" electrical galvanic resistivity or lateral logging 
instrument is shown in Fig. 3. A borehole 30 is shown penetrating an earth formation 
32 in which is suspended a seven-electrode lateral logging instrument 70 by means of 
a wireline cable (not shown) that operates in the same general manner as that described 
in Fig. 2. The survey-current electrode 72 (AO) is disposed in the center of the 
longitudinal array of electrodes, with a pair of measure electrodes 74 (Ml) and 76 
(Ml') being disposed above the center electrode 72 (AO), and a pair of measure 
electrodes 78 (M2) and 80 (M2') being disposed below the center electrode 72 (AO). 
An upper focusing or guard electrode 82 (Al) is disposed above the upper measuring 
electrode 76 (Ml'), while a lower focusing or guard electrode 84 (A2) is disposed 
below the lower measuring electrode 80 (M2'). The pair of upper and lower measure 
electrodes 74 (Ml) and 78 (M2) are interconnected and short-circuited, as are electrode 
pairs 76 (Ml') and 80 (M2') and electrode pairs 82 (Al) and 84 (A2). The guard 
electrodes 82 (Al) and 84 (A2) inject an auxiliary or "bucking" current into the 
formation 32, and the survey electrode 72 (AO) emits a "survey" current of the same 
polarity as the guard electrodes 82 (Al) and 84 (A2). A control signal representing 
the potential difference between the measure electrode pairs 74-78 (M1-M2) and 76-80 
(Ml'-M2') is processed to control the current delivered to the survey electrode 72 
(AO). The system is self-nulling in that the bucking currents produced by the focusing 
electrodes 82 (Al)-84 (A2) is continuously adjusted in order to maintain the difference 
in the voltage differential across the measure electrode pairs 74-78 and 76-80 equal to 
zero [V(74-78) - V(76-80) = 0]. This causes a lateral "focusing" effect in the 
formation surrounding the tool where the current from the survey electrode 72 (AO) is 
forced laterally outwardly to deeply penetrate the formation 32 as shown 
diagrammatically at 86. The potential measurement can be made at any one of the 
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measure electrodes 74, 76, 78 or 80 with respect to a remote electrode spaced a great 
distance from the current field surrounding the tool, typically at the earth surface (not 
shown). The survey current flowing from electrode 72 (AO) is measured, and knowing 
the measured potential then permits calculation of a resistivity measurement R. 

In current field practice, the induction and galvanic (lateral) logging 
measurements are taken separately in different borehole environments (mud 
conductivity) and predicted formation connate water resistivity ratios with the predicted 
resistivity of the mud filtrate in the formation. Generally, when highly conductive 
muds are used (salt water-based mud) the lateral logging instruments are preferred 
since the mud conductivity will provide direct current flow into the formation. On the 
other hand, the induction logging instrument was designed to be used in oil-based or 
low conductive water muds or in air-filled boreholes. Fig. 4 is a standard diagram 
showing the environmental and formation conditions usually considered for determining 
whether to select an induction logging instrument or a lateral logging instrument when 
logging a specific well. For instance, the resistivity of the formation water R,, would 
be estimated using, for example, information from a nearby well. The resistivity of 
the mud filtrate R mf will then be estimated, which can also be done from information 
from a nearby well if the mud system will be the same, or by measurement of a sample 
if one is available. Then the porosity of the formation is estimated, which can be 
based on prior known information of the porosity obtained from a nearby well or from 
other logs that measure porosity. 

Using the estimated data for R w and R^, a calculation of the ratio IWRw is 
made. By referring to a chart or diagram similar to the one shown in Fig. 4, it is 
possible to determine the type of instrument most suited for the particular well. If the 
Rnrf/R^ ratio is 2.5 to 7.5, formation porosity is 15-25% and R w is equal to or greater 
than 1 ohm-m, then the locus of those points will fall above the curve 88 (1^=1 ohm- 
m) and to the right of the vertical dividing line 89, and thus the induction log is 
preferred. Similarly, if the Rmf/R*, ratio is greater than 2.5, formation porosity is 5- 
10% and is equal to or greater than 0.01 ohm-m but less than 0.1 ohm-m, or equal 
to or greater than 0. 1 ohm-m but less than 1 ohm-m, the locus of the data points will 
also fall above curves 88' and 88" and to the right of line 89, again indicating that an 
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induction log is preferred. However, if the R^I^ ratio is below 2.5, and particularly 
if it is below 1.5, the locus of the data points will be to the left of line 89, and below 
all curves 88, 88' and 88", and then the lateral log is preferred. Although both 
instruments may be run in some wells, they are rarely interpreted together because 
either one of the logs is currently deemed to exhibit the most applicable response 
characteristics for each particular formation of interest in the borehole and the 
formation depending on the characteristics referred to above. 

However, the method contemplated herein utilizes the combined responses from 
induction (electromagnetic) and lateral (galvanic) logging instruments to obtain 
improved definition of the electrical parameters of the formation around the borehole. 
The physical principle behind this method can be explained by first looking at the two 
individual principles of resistivity investigation of earth formations, namely "inductive" 
coupling to the formation and measurement of magnetic fields created in the formation, 
and "galvanic" direct electrical coupling to the formation and the measuring of 
electrical fields created thereby. The inductive approach generates induction currents 
in the formation which predominantly flow in conductive formation strata. This means 
that inductive measuring techniques are mainly sensitive to conductive strata in 
formations and relatively insensitive to more resistive strata. On the other hand, the 
galvanic techniques (lateral logging) generate electrical currents in the formation strata 
which in turn generate electrical charges at the layer boundaries of the formation strata. 
If the formation strata medium between two boundaries causing current flow is 
conductive, it will act to short circuit between the resistive layers. Thus the lateral 
logging techniques tend to be more sensitive to resistive strata because the conductive 
layers of the strata (which have been short circuited) do not significantly contribute to 
measuring the formation resistance. This above described measurement phenomena 
also explains why inductive and lateral logging techniques have been commonly applied 
separately to different resistivity (conductivity) formation environments. 

The combination measurement method herein disclosed is based on the 
technique of combining the physical principles of inductive and galvanic measurement 
of resistivity (conductivity) either in one tool or separate tools. The key to the method 
is the combination of the different physical principles of inductive and galvanic 
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measurement in the interpretation of the logging data for forcing theoretical responses 
(including both inductive and galvanic logging responses) of a single model of the 
formation to match real data which will improve the accuracy of determining the 
formation resistivity (conductivity) . 

Such measurements can be derived by using borehole logging devices that 
provide both induction logging and lateral logging types of responses as follows: 

1. two or more separate tools (galvanic and induction) could be used in 
either separate logging runs or as a combined tool string; or 

2. one combined tool measuring the inductive and galvanic response of the 
medium in parallel (simultaneously or separated by a short interval in 
the tool string); or 

3 . separate tools run separately even from different tool families as long as 
the responses characterize the induction and galvanic responses of the 
formation. 

The combination of the responses (measurements) uses the same mathematical 
model for both types of response data sets and the same iterative optimization 
technique(s) which yield one optimized formation model which is consistent with both 
data sets for a given borehole/formation situation. The uncertainty in the model 
parameters is improved through the combination of using two (or more) different data 
sets based on the galvanic and induction responses. The combination may be 
accomplished by using any one of several well known techniques, such as iterative 
numerical methods, geographical techniques, or even direct imaging techniques. The 
key determinant for the selection of the method is the stability and speed of 
convergence to the final formation model. In addition, the reliability of the 
improvements can be evaluated by using error boundary or confidence levels or 
sensitivities to the model parameters. 

To illustrate the preferred embodiment of the method, an example model 
formation will be assumed as shown in Fig. 5. A borehole 90 of diameter D bh is 
shown penetrating a formation of interest 92. The "invaded zone" of the formation 
invaded by the drilling mud fluid is shown at 94, which, in this example, has a step 
profile of diameter Dj. The resistivity of the drilling mud itself in the borehole 90 is 
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shown as R,,,, the resistivity of the invaded zone as and the resistivity of the 
formation as J^. For this example, the borehole diameter D bh was assumed to be 8 
inches and the diameter of the invaded zone (Dj) is 60 inches. All possible cases of 
simple borehole resistivity models can be divided into four (4) basic classes of models 
of resistivity profile: 

1. A-type: in which the resistivity increases from through R xo and to 

2. Q-type: in which the resistivity decreases from through R, 0 and to 

3. K-type: in which the resistivity of the interbedded invaded zone R xo is 
higher than R^ and R t ; and 

4. H-type: in which the resistivity of the interbedded invaded zone R^ is 
lower than R^ and R,. 

The four classes of resistivity profile described above can be summarized in 

the following example: 





Em 


Exo 




A-type 


1 < 


10 < 


100 


Q-type 


10 > 


1 > 


0.1 


K-type 


1 < 


10 > 


1 


H-type 


10 > 


1 < 


100 



Evaluation of all of the above classes in the model would cover to first order all 
possible situations occurring in nature. Using the above established earth model 
parameters and the four classes of resistivity models, predicted responses can be 
calculated for the current effects in the formation due to each of the induction and 
lateral logging types of instruments using Maxwell's equations. 

The next step is to select several known induction and lateral logging tools that 
can be used in various combinations to obtain induction and lateral logging responses 
for each of the four classes of the assumed earth model formations. In the example 
situation herein under explanation, the following types of induction and lateral logging 
instruments of Atlas Wireline Services Division of Western Atlas International, Inc. 
were selected to determine their respective responses to the model formation: 
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1. Dual Phase Induction Logging (DPIL) Instrument (Series 1507), 
including in combination a Shallow Focused Lateral (SFL) Logging 
Instrument; 

2. Dual Laterolog (DLL) Logging Instrument (Series 1229); 

3 . Lateral Logging Instrument having a high vertical resolution and deep 
lateral penetration known as a Thin-Bed Resistivity (TBRT) Logging 

■ 

Instrument (Series 1227); and 

4. Micro Laterolog (MLL) Logging Instrument (Series 1233). (The MIX 
may be run in combination with the TBRT as one tool (the Series 1227) 
as illustrated in Fig. 12.) 

A basic description of the operation and response of each of the above-described 
instruments will follow in connection with Figs. 1-14. 

The basic operation of the DPIL and SFL instruments identified above as used 
in the development of the data for this invention is as described above in connection 
with Figs. 1 and 2, and is well known in the logging art. As above-described, the 
selected combination of the DPIL and SFL instruments may be run in a borehole, 
generating the typical resistivity log as shown in Fig. 6. The DPIL deep response 
curve is shown at 100, the DPIL medium response curve is shown at 102 and the SFL 
shallow focused response curve is shown at 104. The value of R, maybe obtained from 
standard "tornado" charts, well known in the logging art, and as shown in Figs. 7 and 
8. The chart of Fig. 7 is for the situation where R t <R X0 , while the chart of Fig. 8 is 
for the situation where R t >R X0 . The charts of Figs. 7 and 8 be used to compute logs 
based on the following equations: 



Rfl = Jfl x R TO + (1-Jn.) xR t 



(1) 



1/Ru.m =G ILM /R*o + (l-Gn.J/R, 



(2) 



1/R.ld = (G ILD /RJ + (l-G^/R, 



(3) 
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where: R, 0 = resistivity of the formation invaded by the drilling 

mud; 

= resistivity of the undisturbed formation; 

J = pseudo geometric factor for the SFL at the 

invasion diameter; 

G = geometric factor for the DPIL at the invasion 

diameter; 

FL = SFL shallow log; 

ILM = DPIL medium log; and 

ILD = DPIL deep log 
The basic operation of the DLL instrument identified above as used in the 
development of the data for this invention is as described above in connection with 
Figs. 2 and 3, and is well known in the logging art. As above-described, the selected 
DLL instrument may be run in a borehole and will generate the typical resistivity log 
as shown in Fig. 9. The DLL deep lateral response curve is identified at 108, while 
the DLL shallow lateral response curve is shown at 110. For R t >R X0 , the value of ^ 
may be obtained from the example "tornado" chart shown in Fig. 10 in a manner 
similar to the manner in which the value of Rt was obtained for the DPIL tool as 
discussed above and which is well known in the logging art. 

The micro laterolog logging (MLL) instrument identified above operates 
similarly to a shallow lateral logging instrument and is contained in a pad 124 attached 
to an extendable arm 122 from a tool mandrel 116 as shown in Fig. 12. The arm 122, 
when extended, forces the pad 124 into contact with the borehole wall. The lateral 
focusing is created as a result of the individual "button" electrodes 125 and the current 
is forced laterally into the formation, similar to the lateral focusing shown in Fig. 2, 
but directionally from the pad 124 in order to measure a small volume of the formation 
adjacent to the borehole wall and is well known in the logging art. As above- 
described, the selected MLL instrument may be run in a borehole and will generate the 
typical resistivity log as shown in Fig. 11 at 112, which in the example shown, was 
taken simultaneously with a standard DLL instrument. The DLL deep lateral response 
curve is identified at 108*, the DLL shallow lateral response curve is shown at 110', 
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while the MLL response curve is shown at 112. The value of R,,, may be obtained 
from a MLL response in a manner similar to the manner in which the value of R< was 
obtained for the DPIL and DLL instruments as discussed above and well known in the 
logging art. 

Fig. 12 shows a typical tool mandrel 116 for running the TBRT and the MLL 
instruments above described. The mandrel 116 includes a pair of extendable arms 118 
and 122 for carrying pads 120 and 124, respectively. The pad 120 may be the TBRT 
instrument which is a high-resolution lateral logging instrument. The lateral focusing 
is created as a result of using the classic "Laterolog-3" (or Guard Log) design well 
known in the art. The pad 120 design uses a "button" electrode 123 design well known 
in the art by utilizing long "guard" electrodes which focus the survey current into the 
formation, similar to the lateral focusing shown in Fig. 2, but directionally from the 
pad 120 in order to more deeply measure a small volume of the formation adjacent to 
the borehole. Vertical resolution is determined primarily by the survey electrode 
vertical length which can be sufficiently short to achieve the desired vertical 
resolution. The MLL instrument pad 124 operates as hereinabove previously 
described. A typical resistivity log for the TBRT is shown in Fig. 13. The TBRT 
response curve is shown at 130. Fig. 14 is an example of a correction chart for 
borehole corrections to TBRT log readings, where the curves 134 are labelled for 
various borehole diameters. 

Referring to Figs. 15-17, diagrammatic views of different combinations of 
volumes of the earth formation surrounding the borehole that may be common to 
induction and galvanic instruments are shown. In Fig. 15, a borehole 135 is shown 
penetrating an earth formation 132, and a shaded area 136, representing a three- 
dimensional cylindrical volume surrounding the borehole (having a radius v x and a 
height h,), which may be the volume investigated by a selected galvanic instrument. 
The shaded volume 138-138' radially displaced from the borehole, represents a three- 
dimensional "doughnut" shaped volume surrounding the borehole (having an outer 
radius of r 2 and a height dimension hj), may be the volume investigated by a selected 
induction instrument. The two volumes 136 and 138-138' overlap in a common 
volume 140 (represented by the cross-hatched area) that is investigated by both the 
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galvanic and induction instruments. 

Fig. 16 shows the earth formation 132 penetrated by a borehole 135 surrounded 
by a radially extending shaded area 137, representing a three-dimensional cylindrical 
volume surrounding the borehole (having a radius of r 2 and a height hj) which may 
represent a volume investigated by a selected galvanic instrument having a greater 
radius of penetration than the galvanic instrument producing the volume 136 shown in 
Fig. 15. The shaded "doughnut-shaped volume 139-139' (having an outer radius of 
r 2 and a height dimension of h 2 ) represents one radially spaced three-dimensional 
volume surrounding the borehole which may be investigated by a selected induction 
instrument. Similarly, the shaded "doughnut"-shaped volume 141-141* (having an 
outer radius of r 3 and a height of h 2 , approximately the height of the volume 139-139') 
represents a second radially spaced three-dimensional volume surrounding the borehole 
and radially spaced from the first such volume 139-139' and may also be the volume 
investigated by a second selected induction instrument. The volumes 137, 139-139' 
and 141-141' form overlapping common volumes 142 and 144 represented by the 
"cross-hatched" areas having a height hj and the outer radii of r, and r 2 , respectively. 
The volumes 139-139' and 141-141' are radially displaced from the borehole at 
different radial distances and constitute common volumes of formation investigated by 
the combination of galvanic and induction logging instruments. 

In Fig. 17, the earth formation 132 is shown penetrated by a borehole 135 and 
surrounded by a pair of radially extending and vertically-spaced shaded areas 143 and 
145 that each represents a three-dimensional cylindrical volume surrounding the 

9 

borehole and having a radius x height combination r 2 x h 2 and r, x h u respectively, 
where h, and h 2 are vertically spaced apart along the borehole, and which may 
represent a volume investigated by a pair of different selected galvanic instruments or 
a single galvanic instrument having differing depths of penetration. The shaded 
"doughnut "-shaped volume 146-146' (having an outer radius of r 3 and a height of h 3 ) 
represents one radially spaced three-dimensional volume surrounding the borehole 
which may be investigated by a selected induction instrument. The different volumes 
143, 145 and 146-146' form common overlapping volumes 147-147' and 149-149', 
represented by the "cross-hatched" areas having an outer radius r 3 and spaced-apart 
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vertical heights h 2 and hj respectively. The volumes 147-147' and 149-149' are 
radially spaced about the borehole and vertically spaced apart at different distances and 
constitute a common volume of formation investigated by the combination of galvanic 
and induction logging instruments. 

The next step in this illustration of the preferred embodiment of the method of 
combining galvanic and induction tool responses is the depiction of the predicted 
responses of combinations of the above described induction and lateral logging 
instruments, which is accomplished by calculating ellipsoids of confidence that form 
a "three-dimensional" surface of the type shown at 150 in Fig. 18. The technique of 
calculating the ellipsoids of confidence as shown at 150 will now be described. 

If all of the n measurements are put in a sequence as follows: 



f= f, 



(4) 



and all of the m model parameters are put in a sequence as follows: 
P = Pi 

(5) 

Pa. 

then the components of vector f are functions of vector jTas follows: 

? = T (p) ( 6) 

If p 0 is an estimated (a priori) model and f 0 is a response from this model, 
then the difference ( f - f 0 ) may be represented in the form: 

A. ^ 

(f-fj = D (p- p 0 ) (7) 



where D is a matrix n x m of partial derivatives of f with respect to 
p as follows: 

df; 



Dy = dpj (8) 
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where i = l,n; 

j=l,m; 

The same formula for (f-f 0 ) can be represented in terms of relative 
deviations: 

5f = D 6p 

f -f 



where: rj 

6t { = 



foi 



(9) 



Spj = PfPoj (10) 
Pj 

Dij = 3£nfi 



If data errors are distributed normally by calculating X 2 - statistics of 6f, we 
obtain the following: 

X 2 (?f T f 5f)= X?(5 r 'lP T fD] 6p) (12) 

where: E is the matrix of the second moments of the random vector df . 
It follows from equation (12) that the surfaces of constant values of X 2 statistics 
argument are the surfaces under consideration which are ellipsoidal in shape: 

(5p T D T f D 6p T ) = [X 2 ] 1 (c) (13) 
where c is a given confidence level. 

For the modelling herein described, the assumption is made that the noise is 
randomly distributed and the data is measured with an accuracy of one percent (1,0%). 
It is also required that the confidence level (c) of the interpreted model parameter is 
at ninety-five percent (95.0%). The axes of the three-dimensional coordinate system 
shown in Fig. 18 are labeled for the percentage deviations from the earth model 
parameters values R xo , Rt and D s . The mud resistivity is assumed to be known since 
it is routinely measured and corrected for temperature when applying standard logging 
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techniques. As shown in the example of Fig. 18, the ellipsoids of confidence really 
form a three-dimensional surface 150, and it is hard to visualize the range of the 
confidence levels or values. Accordingly, if the ellipsoids of confidence of the surface 
150 are projected onto each of the three axes of the percentage deviations of the 
parameters R x0 , Rt and D { in Fig. 18, the projected images will be as seen as projected 
two-dimensional ellipsoid images 152 along axis (AD/Dj, AR X0 /R, 0 ), 154 along axis 
(AR/R,, AR^/Rjo) and 156 (AD/Dj, AltyR,). From such projections, one can visualize 
the range of uncertainties in the interpreted parameters. 

Assuming the A-type model of formation as previously described, with the 
following known parameters: R m = 1, R^ = 10, R t = 100, D bh = 8 in., and D ; = 60 
in. , it has been found useful that comparison error ellipsoids may be determined for 
the following combinations of the selected induction and lateral logging instruments 
hereinabove described: 

1. Combination of DPIL + SFL; 

2. Combination of DPIL + SFL + MLL + TBRT; 

3. Combination of DLL + TBRT + MLL; and 

4. Combination of DPIL + SFL + DLL + MLL + TBRT. 

Of course, other combinations may be used, and the invention is certainly not restricted 
to the above listed or described combinations of induction and lateral logging 
instruments and responses. However, the above-described combinations were used and 
such combinations achieved the desired results. 

Using a joint inversion mathematical interpretation technique, the error ellipsoids 
for the A-type model formation and the four (4) combinations of induction and lateral 
logging instruments as above described are calculated. The "joint inversion" technique 
used is preferably of the type disclosed by Strack, "Exploration with Deep Transient 
Electromagnetics," Elsevier, 1993, which technique is hereby incorporated for all 
purposes. The comparisons of the error ellipsoids for the A-type model formation and 
the four tool combinations are shown in Fig. 19. The combination of the DPIL and 
the SFL responses have the largest error as shown by the ellipsoid projections for the 
combination of the DPIL (induction response) + the SFL (shallow galvanic logging 
response) shown as curves 160, 162 and 164. Combining these responses with the 
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responses of the TBRT and MLL (galvanic logging instruments) significantly improves 
the error situation, as shown by curves 166, 168 and 170. However, if the induction 
instrument (DPIL) and spherically focused lateral logging responses (SFL) are removed 

» 

and the dual lateral logging (DLL) responses are added, the error gets worse as is 
shown by the larger curves 172, 174 and 176. On the other hand, if the responses of 
all of the instruments are combined, the error curves are substantially reduced as 
shown by the curves 178, 180 and 182 (the curve 182 is substantially covered by and 
obstructed by the larger curve 176). The improvement of error varies with the 
measurement of the various parameters by factors of 3 to 15. 

The information shown in the diagram of Fig. 19 can also be displayed in the 
form of a logarithmic model plot as shown in Fig. 20, where the true model is depicted 
by the single line towards the center of the suite of curves for each of the induction 
and/or lateral logging tool combinations. The measurement uncertainties as resulting 
from each of the various combinations are represented by pairs of lines either side of 
the true value. The combination for the DPIL + SFL is shown by the uncertainty 
envelopes defined by the curves 186 and 186'. The uncertainty envelope for the 
combination of the DPIL + SFL + MLL is shown by the curves 188 and 188'. The 
uncertainty envelope for the combination of the DLL + TBRT + MLL is shown by 
the curves 190 and 190' , while the uncertainty envelope for the combination of all tools 
is shown by the curves 192 and 192'. 

If it is assumed that the Q-type model formation has the following known 
parameters: ^ = 10, R ro = 1, R, = 0.1, D bh = 8 in. and D t = 60 in., then using the 
four (4) combinations of induction and galvanic logging instruments as hereinabove 
described, the joint inversion technique will calculate the error ellipsoids for each such 
combination in the same manner as above described for the A-type model formation. 
The comparisons of the error ellipsoids for the Q-type model formation and the four 
tool combinations are shown in Fig. 21. The combination of the DPIL + SFL 
responses (primarily the inductive response) have a large error as reflected by the 
ellipsoid projections for the combination of the DPIL (induction response) + the SFL 
(shallow galvanic logging response) shown as curves 194, 196 and 198. Again, 
combining these responses with the responses of the TBRT + MLL (galvanic logging 
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instruments) significantly improves the error situation, as shown by curves 200, 202 
and 204. However, if the DPIL + SFL responses are removed and the dual lateral 
logging (DLL) responses are added (substantially only lateral focused responses), the 
error is even greater as is shown by the larger curves 206, 208 and 210. As is true 
for the A-type model formation, if the responses of all of the instruments are 
combined, the error curves are substantially reduced as shown by the curves 212, 214 
and 216. 

The information shown in the diagram of Fig. 21 can also be displayed in the 
form of a logarithmic model plot as shown in Fig. 22, where the true model is always 
towards the center of the suite of curves for each of the induction and/or lateral logging 
tool combinations. The combination for the DPIL + SFL is shown by the uncertainty 
envelopes defined by the curves 220 and 220'. The uncertainty envelope for the 
combination of the DPIL + SFL + MLL + TBRT is shown by the curves 222 and 
222'. The uncertainty envelope for the combination of the DLL + TBRT + MLL is 
shown by the curves 224 and 224' , while the uncertainty envelope for the combination 
of all tools is shown by the curves 226 and 226*. 

In the third case, assume the K-type model formation has the following known 
parameters: 1^=1, R xo = 10, = 1, D bh = 8 in. and Dj = 60 in., and the four (4) 
combinations of induction and galvanic logging instruments are used as hereinabove 
described, then the error ellipsoids for each such combination can be calculated in the 
same manner as above described for the A-type model formation. The comparisons 
of the error ellipsoids for the K-type model formation and the four tool combinations 
are shown in Fig. 23. The combination of the DPIL + SFL responses (primarily the 
inductive response) has the largest error as reflected by the ellipsoid projections for the 
combination of the DPIL (induction response) + the SFL (shallow galvanic logging 
response) shown as curves 228, 230 and 232. If the first combination is further 
combined with the responses of the MLL + TBRT (lateral logging instruments) the 
error situation is significantly improved, as shown by curves 234, 236 and 238. 
However, if the DPIL + SFL responses are removed and the dual lateral logging 
(DLL) responses added (substantially only lateral focused responses), the error is even 
greater as is shown by the larger curves 240, 242 and 244. If the responses of all of 
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the instruments are combined, the error curves are substantially reduced as shown by 
the curves 246, 248 and 250, which are substantially coincident with the curves 234, 
236 and 238 of the combination of the DPIL + SFL + MLL + TBRT instrument 
responses. 

As previously shown in Figs. 20 and 21, the information shown in the diagram 
of Fig. 23 for the K-type of formation model can also be displayed in the form of a 
logarithmic model plot as shown in Fig. 21, where the true model is always towards 
the center of the suite of curves for each of the induction and/or galvanic logging tool 
combinations. The combination for the DPIL + SFL is shown by the largest 
uncertainty envelope defined by the curves 254 and 254'. The uncertainty envelope 
for the combination of the DPIL + SFL + MLL + TBRT is shown by the curves 256 
and 256\ The uncertainty envelope for the combination of the DLL + TBRT + MLL 
is shown by the curves 258 and 258', which shows a much enlarged error of 
uncertainty. The uncertainty envelope of the combination of all tools is shown by the 
curves 260 and 260', which are seen to virtually coincide with the curves 256 and 256* 
of the combination of the DPIL + SFL + MLL + TBRT. 

In the last case, assuming that the H-type model formation has the following 
known parameters: R™ = 10, R x0 = 1, R, = 100, D bh = 8 in. and D { = 60 in., and 
that the four (4) combinations of induction and lateral logging instruments are used as 
hereinabove described, then the error ellipsoids for each such combination can be 
calculated in the same manner as above described for the other type of model 
formations . The comparisons of the error ellipsoids for the H-type model formation 
and the four tool combinations are shown in Fig. 25. The combination of the DPIL 
+ SFL responses (primarily the inductive response) has the largest error as reflected 
by the ellipsoid projections for the combination of the DPIL (induction response) + the 
SFL (shallow galvanic logging response) shown as curves 264, 266 and 268. Further 
combining the first combination with the responses of the MLL + TBRT (galvanic 
logging instruments) the error situation is significantly improved, as shown by curves 
270 (the curve 270 is substantially covered by the larger curves 276 and 282 for the 
other combinations as will hereinafter be described), 272 and 274. As in the earlier 
cases, if the DPIL + SFL responses are removed and the dual lateral logging (DLL) 
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responses added (substantially only lateral focused responses), the error is greater as 
is shown by the larger curves 276, 278 and 280. Combining all of the responses for 
all of the instruments causes all of the error curves to be substantially reduced as 
shown by the curves 282, 284 and 286. The curves 284 and 286 are substantially 
coincident with the larger curves 278 and 280 of the combination of the DLL + TBRT 
+ MLL and are obscured in Fig. 25. 

As previously shown in Figs. 20, 22 and 24, the information shown in the 
diagram of Fig. 25 can also be displayed in the form of a logarithmic model plot as 
shown in Fig. 26, where the true model is always towards the center of the suite of 
curves for each of the induction and/or lateral logging tool combinations. The 
combination for the DPIL + SFL is shown by the largest uncertainty envelope defined 
by the curves 290 and 290' . The uncertainty envelope for the combination of the DPIL 
+ SFL + MLL + TBRT is shown by the curves 292 and 292' . The uncertainty 
envelope for the combination of the DLL + TBRT + MLL is shown by the curves 
294 and 294', while the uncertainty envelope of the combination of all tools is shown 
by the curves 296 and 296', which again can be seen to virtually coincide with the 

curves 294 and 294*. 

In practice, the combination of the two principles of physical measurement 
(inductive and galvanic) is accomplished by combining the independent data sets from 
each measurement and deriving a single model which gives responses in agreement 
with the data sets from each type of measurement. The determination of the agreement 
of the data sets with the derived model can be expressed, for example, as an RMS 
error or fit. The interpretative technique for deriving this model is well known as 
"inversion" or "least square inversion" or "system identification" in the data processing 
field. Also, neural network techniques or other artificial intelligence tools could be 
used to achieve the same result. The processing technique used in initially deriving the 
invention as above described uses the derivatives of the model function with respect to 
the modeled parameters to evaluate the sensitivities of the modelled solution to each 
of the parameters. This approach facilitates the determination of which parameters are 
sensitive due to the induction logging and lateral logging measurements, respectively. 
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In the field application of this technique, the "joint inversion" data processing 
technique, as hereinabove described, will be used and a functional block diagram of 
an example process 300 is shown in Fig. 24. As previously described, the known 
available ("a priori") geological, physical and environmental information is also 
supplied at functional step 302, which, combined with the selected induction and 
galvanic tool parameter data, is then used to construct an initial assumed formation 
model at functional step 304, similar to one of the example type of formations (A-, Q-, 
K- or H-type) or combinations thereof, as above described. This initial model, 
determined at step 304, becomes the earth model shown at step 306. The information 
to do this may be compiled from multiple and complex origins, or simply be an 
analysis of the data from one of the induction or galvanic devices providing the data 
set. Based on the models for the selected induction and galvanic tools as determined 
at step 308, the estimated parameters of the initial earth model as determined at the 
functional step 306 are then used to calculate the selected inductive and lateral logging 
instrument response functions, as hereinabove described and, as shown by the dotted 
lines 310 and to generate "synthetic" modelled data at step 312 for the initial earth 
model determined at step 306. 

The selected inductive and galvanic logging instruments are used in the field 
as shown in step 314 to interact with the actual earth formation surrounding the 
borehole as shown at 316 to obtain actual field data at step 320. The field data 
obtained in functional step 320 and the synthetic modelled data obtained in functional 
step 312 are compared at step 322 and the deviation is applied at step 324 to compare 
the deviation with an established range of values of resemblance criteria between the 
field measurement response data and the calculated measurement modelled response 
data. During the first iteration, the comparison of the actual field data obtained by the 
selected inductive and electrical galvanic logging instruments at step 320 to the 
"synthetic" or "modelled" data generated by the initial earth model and the selected 
logging instrument responses as obtained at step 312 is likely to differ by an amount 
greater than an established range of values of a resemblance criteria between the field 
measurement response data and the calculated measurement data. Any deviation values 
or data that occurs at the comparison step 324 is used at the functional step 326 to 
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provide adjustment data for the initial earth model formation parameters via feedback 
loop 327 to the earth formation model step 306 using standard data processing 
techniques. 

The adjusted earth formation model 306 developed by the model change step 
326 is again used to generate "synthetic" or "modelled" instrument response data based 
on the selected inductive and galvanic logging instruments as shown at 308. Such 
adjusted "synthetic" or "modelled" logging response data generated at step 312 is again 
compared with the actual combined field logging data obtained in step 320, using the 
same selected inductive and electrical galvanic logging instruments, as hereinabove 
described. If the numerical value of the resemblance criteria is too large at the 
comparison step 324 in this second comparison or iteration, the earth formation model 
is again adjusted at step 326, as hereinabove described, and the iterative process 
continues by the repetition of functional steps 312, 322, 324 and 326. When the 
comparison made at step 322 finally derives a deviation that falls within the 
predetermined criteria of resemblance at step 324, the adjusted formation model is 
considered to be optimized as shown at step 328. Based on the optimization of the 
adjusted model at step 328, the actual combined field data obtained at step 320 can now 
be applied as shown at 330 to the final model obtained in step 328 to generate the 
desired "optimized" output parameter, as shown at functional step 332. If the desired 
formation output parameter is a "resistivity" measurement of the formation strata of 
interest, then the output parameter of functional step 332 will typically include measure 
of R,. However, if the desired formation output parameter is a "conductivity" 
measurement of the formation strata of interest, then the output parameter of functional 
step 332 will typically include a measure of C t . The "inversion" modelling technique 
334 previously referred to encompasses the functional steps included within the dashed 
lineblock and includes steps 306, 308, 310, 312, 322, 324, 326 and 328. 

The examples described above demonstrate that the technique of combining the 
induction and galvanic effects of induction and lateral logging instruments is at least 
as good as and generally superior to and more accurate than relying on the 
interpretation of responses to any single tool. The method reduces the ambiguity in 
inteipretation for a particular formation using any individual resistivity logging 
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instrument. While it is possible that there may exist some instances where the 
formation model parameters might not permit this method to determine the optimum 
resistivity for the formation in question, the use of this optimization method when 
combined with other interpretative techniques as well as additional geological data will 
help to obtain better final results. In addition, when applying this technique to more 
complicated situations and models, certain combinations of the model parameters could 
be the final resolved output of the method as determined at step 328 rather than the 
independent output obtained at step 332. 

It is common as part of the well logging art to measure the porosity of the 
formation, 0 , and to calculate the water saturation (S w ) of the formation of interest 
using a variation of the well known Archie equation: 

where: S w = water saturation of the formation 

= resistivity of the connate water of the formation 
1^ = resistivity of the formation of interest 
0 = formation porosity 

c = constant; (1.0 for carbonates & 0.90 for sands) 

Once the water saturation, S w , has been determined, the hydrocarbon saturation, 
S h , of the formation can readily be calculated by 

S h = (1-S W ) (15) 
The fraction of total formation volume that is hydrocarbons can be calculated as 0(S h ) 
or 

0 (1-SJ. 

As may readily be seen and appreciated, using the herein described method to 
more accurately measure formation parameters, especially the resistivity of the 
formation, R t , using the combination of galvanic and induction methods, will provide 
more accurate determination of the formation hydrocarbon saturation, S h , and 
ultimately, the key calculation, the quantity of hydrocarbons present in the formation. 

Numerous variations and modifications may be made in the structure and/or 



WO 95/03557 



PCT/US94/07991 



-27- 

steps of the method and/or process herein described without departing from the present 
invention. Accordingly, it should be clearly understood that the forms of the invention 
herein described and shown in the figures of the accompanying drawings are illustrative 
only and are not intended to limit the scope of the invention. Practical application of 
the invention might include simplified processing to allow quick estimates of results to 
be made in real time or near real time. 
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CLAIMS 

1. A method of determining selected parameters of an earth formation 
surrounding a borehole, comprising the following steps 

(a) obtaining at least one induction logging measurement of the 
selected parameters in a first predetermined volume of the formation surrounding the 
borehole having known first radial and vertical dimensions, 

(b) obtaining at least one galvanic logging measurement of the 
identical selected parameters in a second predetermined volume of the formation 
surrounding the borehole having known second radial and vertical dimensions that 
overlap said first radial and vertical dimensions of said first predetermined volume, 
whereby said overlapping volumes form a representative common volume of the 
formation, and 

(c) combining said induction and galvanic logging measurements to 
obtain a measurement of the selected parameters of the earth formation surrounding the 
borehole in said representative common volume of the formation. 

2. The method of claim 1, wherein said combining step comprises the 
following 

(a) combining the field measurement response data for each of said 
selected inductive and electrical galvanic logging instruments for said representative 
common volume of the formation; 

(b) determining an initial earth formation model based on known and 
estimated geological, physical and environmental information; 

(c) calculating for said initial earth formation model and the known 
selected inductive and electrical galvanic logging instruments the anticipated formation 
responses; 
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(d) determining from said formation responses the deviation between 
the formation modelled data and the measured data; 

(e) establishing a range of values for the criteria of resemblance 
between said field measurement response data and said calculated measurement data; 

(f) adjusting the initial earth formation model in response to said 
deviation between said anticipated formation response data and said combined field 
measurement data exceeding said established criteria of resemblance determined in step 

(e); 

(g) repeating steps (c), (d) and (0 until said deviation between said 
calculated earth formation response and said actual combined field data is within said 
established criteria of resemblance as determined in step (e); and 

(h) in response to step (g) calculating an output value representative 

* 

of the measurement of the selected earth formation parameters. 

* 

3. A method of determining selected parameters of an earth formation 
surrounding a borehole, comprising the following steps 

(a) obtaining at least one induction logging measurement of the 
selected parameters in a first predetermined volume of the formation surrounding the 
borehole, such first predetermined volume having known first radial and vertical 
dimensions, 

(b) obtaining at least one galvanic logging measurement of the 
identical selected parameters in a second predetermined volume of the formation 
surrounding the borehole, such second predetermined volume having known second 
radial and vertical dimensions that differ from said first radial and vertical dimensions 
of said first predetermined volume, whereby said first and second predetermined 
volumes having said known first and second radial and vertical dimensions form a 
representative sample of the surrounding formation that includes multiple common 
volumes having common overlapping radial and vertical dimensions and which are 
radially or vertically spaced from each other, 

(c) combining said induction and galvanic logging measurements to 
obtain a measurement of the selected parameters of the earth formation surrounding the 
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borehole in said representative sample of the surrounding formation that includes 
multiple common volumes having common overlapping radial and vertical dimensions 
and which are radially or vertically spaced from each other. 

4. The method of claim 3, wherein said combining step comprises the 
following 

(a) combining the field measurement response data for each of said 
selected inductive and electrical galvanic logging instruments for said representative 
sample of the surrounding formation; 

(b) determining an initial earth formation model in said representative 
sample of the surrounding formation based on known and estimated geological, 
physical and environmental information; 

(c) calculating for said initial earth formation model and the known 
selected inductive and electrical galvanic logging instruments the anticipated formation 
responses; 

(d) determining from said formation responses the deviation between 
the formation modelled data and the measured data; 

(e) establishing a range of values for the criteria of resemblance 
between said field measurement response data and said calculated measurement data; 

(f) adjusting the initial earth formation model in response to said 
deviation between said anticipated formation response data and said combined field 
measurement data exceeding said established criteria of resemblance determined in step 

(e); 

(g) repeating steps (c), (d) and (f) until said deviation between said 
calculated earth formation response and said actual combined field data is within said 
established criteria of resemblance as determined in step (e); and 

(h) in response to step (g) calculating an output value representative 
of the measurement of the selected earth formation parameters, 

5. A method of determining selected parameters of an earth formation 
surrounding a borehole in accordance with the following steps: 
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(a) selecting inductive logging and electrical galvanic logging 
instruments for measuring responses of the earth formation surrounding the borehole; 

(b) acquiring field measurement response data for each of said 
selected inductive and electrical galvanic logging instruments for the earth formation 
to be measured; 

(c) combining the field measurement response data for each of said 
selected inductive and electrical galvanic logging instruments; 

(d) determining an initial earth formation model based on known and 
estimated geological, physical and environmental information; 

(e) calculating for said initial earth formation model and the known 
selected inductive and electrical galvanic logging instruments the anticipated formation 
responses; 

(f) determining from said formation responses the deviation between the 
formation modelled data and the measured data; 

(g) establishing a range of values for the criteria of resemblance 
between said field measurement response data and said calculated measurement data; 

(h) adjusting the initial earth formation model in response to said 
deviation between said anticipated formation response data and said combined field 
measurement data exceeding said established criteria of resemblance determined in step 

(g); 

(i) repeating steps (e), (f) and (h) until said deviation between said 
calculated earth formation response and said actual combined field data is within said 
established criteria of resemblance as determined in step (g); and 

(j) in response to step (i) calculating an output value representative of 
the measurement of the selected earth formation parameters. 

6. The method as described in claim 5, wherein the earth formation 
surrounding a borehole has a selected volume determined by the following steps 

(a) obtaining at least one induction logging measurement of the 
selected parameters in a first predetermined volume of the formation surrounding the 
borehole having known first radial and vertical dimensions, and 
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(b) obtaining at least one galvanic logging measurement of the 
identical selected parameters in a second predetermined volume of the formation 
surrounding the borehole having known second radial and vertical dimensions that 
overlap said first radial and vertical dimensions of said first predetermined volume, 
whereby said overlapping volumes form a selected common volume of the formation. 

7. The method as described in claim 5, wherein the earth formation 
surrounding a borehole has a selected volume determined by the following steps 

(a) obtaining at least one induction logging measurement of the 
selected parameters in a first predetermined volume of the formation surrounding the 
borehole, such first predetermined volume having known first radial and vertical 
dimensions, and 

(b) obtaining at least one galvanic logging measurement of the 
identical selected parameters in a second predetermined volume of the formation 
surrounding the borehole, such second predetermined volume having known second 
radial and vertical dimensions that differ from said first radial and vertical dimensions 
of said first predetermined volume, 

whereby said first and second predetermined volumes having said known 
first and second radial and vertical dimensions form a representative sample of the 
surrounding formation that includes multiple common volumes having common 
overlapping radial and vertical dimensions and which are radially or vertically spaced 
from each other. 

8. The method of claims 1 or 3 or 5 wherein said value representative of 
the formation parameter using said combined induction and galvanic measurements of 
the parameter may be obtained with greater accuracy than may be derived from using 
either said induction or galvanic measurements alone. 

9. The method of claims 1 or 3 or 5 wherein said value representative of 
the formation parameter using said combined induction and galvanic measurements of 
the parameter may be used to derive a more accurate calculation of hydrocarbon 
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saturation in the earth formation than may be derived from using either said inductive 
or galvanic measurements alone. 

10. The method of claims 1 or 3 or 5 wherein the selected parameter of the 
formation to be determined is resistivity. 

11. A method of determining resistivity of a preselected volume of earth 
formation surrounding a borehole in accordance with the following steps: 

(a) selecting inductive logging and electrical galvanic logging 
instruments for measuring responses related to resistivity of the earth formation 

surrounding the borehole; 

(b) acquiring field measurement response data related to resistivity for 
each of said selected inductive and electrical galvanic logging instruments for the 
selected volume of earth formation to be measured; 

(c) combining the field measurement response data for each of said 
selected inductive and electrical galvanic logging instruments; 

(d) determining an initial earth formation model for the earth 
formation volume to be measured based on known and estimated geological, physical 
and environmental information; 

(e) calculating for said initial volume of earth formation model and 
the known selected inductive and electrical galvanic logging instruments the anticipated 
formation responses; 

(f) determining from said formation responses the deviation between the 
formation modelled data and the measured data; 

(g) establishing a range of values for the criteria of resemblance 
between said field measurement response data and said calculated measurement data; 

(h) adjusting the initial earth formation model in response to said 
deviation between said anticipated formation response data and said combined field 
measurement data exceeding said established criteria of resemblance determined in step 

(g); 
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(i) repeating steps (e), (f) and (h) until said deviation between said 
calculated earth formation response and said actual combined field data is within said 
established criteria of resemblance as determined in step (g); and 

(j) in response to step (i) calculating an output value representative of 
the measurement of resistivity for the selected volume of earth formation. 

12. The method of determining resistivity as described in claim 11, wherein the 
earth formation surrounding a borehole has a selected volume determined by the 
following steps 

(a) obtaining at least one induction logging measurement of the 
selected parameters in a first predetermined volume of the formation surrounding the 
borehole having known first radial and vertical dimensions, and 

(b) obtaining at least one galvanic logging measurement of the 
identical selected parameters in a second predetermined volume of the formation 
surrounding the borehole having known second radial and vertical dimensions that 
overlap said first radial and vertical dimensions of said first predetermined volume, 
whereby said overlapping volumes form a selected common volume of the formation. 

13. The method of determining resistivity as described in claim 11, wherein the 
earth formation surrounding a borehole has a selected volume determined by the 
following steps 

(a) obtaining at least one induction logging measurement of the 
selected parameters in a first predetermined volume of the formation surrounding the 
borehole, such first predetermined volume having known first radial and vertical 
dimensions, and 

(b) obtaining at least one galvanic logging measurement of the 
identical selected parameters in a second predetermined volume of the formation 
surrounding the borehole, such second predetermined volume having known second 
radial and vertical dimensions that differ from said first radial and vertical dimensions 
of said first predetermined volume, 
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whereby said first and second predetermined volumes having said 
known first and second radial and vertical dimensions form a representative sample of 
the surrounding formation that includes multiple common volumes having common 
overlapping radial and vertical dimensions and which are radially or vertically spaced 
from each other. 
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